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Abstract 

The thermal conductivity of YbRh 2 Si 2 has been measured down to very low temperatures under 
field in the basal plane. An additional channel for heat transport appears below 30 mK, both 
in the antiferromagnetic and paramagnetic states, respectively below and above the critical field 
suppressing the magnetic order. This excludes antiferromagnetic magnons as the origin of this 
additional contribution to thermal conductivity. Moreover, this low temperature contribution 
prevails a definite conclusion on the validity or violation of the Wiedemann-Franz law at the 
field-induced quantum critical point. At high temperature in the paramagnetic state, the thermal 
conductivity is sensitive to ferromagnetic fluctuations, previously observed by NMR or neutron 
scattering and required for the occurrence of the sharp electronic spin resonance fracture. 

PACS numbers: 71.27.+a, 72.15.Eb, 74.40.Kb, 75.50.Ee 
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Traditionally, heavy fermion intermetallics are viewed as systems lying close to the insta¬ 
bility to an ordered magnetic state. The instability is driven by the competition between 
Kondo screening of the local moments, and long range order of these moments, induced 
by RKKY interactions. Along the past decades, experimental and theoretical studies have 
endeavored to unveil the evolution of heavy fermions through the instability, focusing on the 
role of quantum fluctuations around the quantum critical point (QCP) pQ. Presently, the 
most prominent scenarios involve instabilities approaching the critical point from the para¬ 
magnetic (PM) side, either on some “hot spots”, which is the so-called “spin density wave” 
scenario, or extended over the whole surface, for the “local quantum critical scenario”. In 
the last case, the Kondo screening mechanism would breakdown near the quantum critical 
point, leading to a change from a large Fermi surface (including the local moments) in the 
paramagnetic state, to a “small Fermi surface” when Kondo coupling disappears |2j|. As a 
consequence, it is also believed that in this last scenario, which involves a collapse of the 
mechanism driving the heavy quasi-particles formation, when approaching the well-defined 
QCP, the low energy excitations would not be described by quasiparticles anymore. So ex¬ 
periments probing this local quantum critical scenario, are mainly trying to reveal a change 
of Fermi surface (FS) volume across the QCP or a clear violation of predictions based on 
the notion of quasiparticles. 

The antiferromagnetic (AF) heavy fermion YbRh 2 Si 2 , has been a main playground for 
both experimental approaches, and strongly pushed forward as an archetype of heavy fermion 
driven by the local quantum criticality scenario Pi- The proximity of YbRh 2 Si 2 to AF 
quantum criticality is attested by the low value of its Neel temperature Tv ~ 70 mK, far 
smaller than its Kondo temperature Tk ~ 25 K |5|. Experimentally, a major appeal of this 
system is that fine field tuning of the AF-PM border can be used to monitor the proximity 
to the QCP [3j (5] . Due to the large anisotropy of the magnetic properties, the critical field 
H c driving the AF-PM instability is also anisotropic, respectively 0.066 T and 0.6 T in the 
easy magnetization plane (a,b), and along the c axis of the body centered tetragonal crystal. 
Evidence of FS reconstruction at H c are mainly based on Hall effect measurements PP- 
However, heavy fermions are multiband systems, so that the interpretation of the Hall effect 
is not straightforward. Other concomitant proof of the FS reconstruction in YbRh 2 Si 2 must 
be found. No convincing changes of the FS at H c appear on the thermoelectric power Pi- 

Recently, the debate has been focused on the detection of a breakdown of the Fermi-liquid 
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regime at H c , using thermal transport at low temperature as a probe: A violation of the 
Wiedemann-Franz law (WFL) would point out a deep change of the electronic excitations 
responsible for charge and heat transport 0 EU H21- In metallic systems, at T —> 0 K, 
when quasiparticle scattering is dominated by elastic mechanisms, all reported cases but 
one satisfy the WFL, stating that L —y Lq in the low temperature limit, where L = trp/T is 
the Lorenz ratio and Lq = L(T —> 0) is the Lorenz number (Lq = (re 2 /3)(kB/e) 2 )- The only 
case reported so far that would not satisfy the WFL in this limit is CeCoIn 5 , for the heat 
current along the c-axis [I3J, which has been put forward as a demonstration of breakdown 
of the quasiparticle picture in a “hot spot” QCP scenario. But even in this last case, one 
cannot exclude that measurements below 30 mK (the lowest temperature reached in [13]) 
would invalidate the violation of the WFL. 

Three different groups have recently published thermal conductivity (k) experiments 
coupled with careful resistivity (p) experiments in YbRh 2 Si 2 , in order to study the temper¬ 
ature dependence of the Lorenz ratio L(T ) = np/T. Contradictory conclusions have been 
drawn by the three groups. The first one [5] claims that the WFL is violated at H c , as 
L(T —* 0) ~ 0.90Lo, while the others concluded conservatively that L(T —* 0) at H c is L 0 
PHH 2\- However, all the experimental data of the three groups are in good agreement: 
The main differences lie in the interpretation, and in the lowest temperature reached, i.e. 
in the extrapolation of the zero temperature values. This last point is a very sensitive is¬ 
sue, as an additional contribution to the electronic thermal conductivity appears at very 
low temperatures (below 30mK) in YbRh 2 Si 2 , which is very difficult to separate from the 
quasiparticles contribution. The measurements we report in the present paper, were realized 
down to 10 mK, twice lower than in those reported in [5], and with a strong emphasis on 
the analysis of this extra contribution, notably to detect if anything specific happens at the 
critical field H c . 

The single crystal investigated in the present study was grown out of In flux. It is the 
same as used in Ref. ra for the thermoelectric power and magnetoresistance studies. The 
residual resistivity ratio RRR = p(300K)/p(0K) of the crystal is 65, the observation of 
quantum oscillations in the resistivity indicates its high quality. The magnetic field was 
applied parallel to the current (longitudinal configuration: H//jq/ /[1,1,0]), up to 4T. The 
zero field measurements have been performed with a zero field cooled magnet. At low field, 
a calibrated Hall sensor has been used to determine precisely the applied magnetic field. 
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FIG. 1. (color online) Electrical resistivity (full squares) and thermal resistivity (open circles) of 
YbRli 2 Si 2 at OT (in black) and at 70 mT, close to H c (in red) at low temperature. Clear anomalies 
are observed in zero field at Tjy, and a new channel for heat transport appears below 30mK (see 
text). Inset: Normalized ratio L(T)/Lq in the whole temperature range at OT and 70mT. 

The thermal conductivity was measured using a standard two thermometers-one heater 
setup, and was checked to be independent of the heat gradient applied, with A T/T in 
the range 0.5-10 %. Calibration of the thermometers was realized in situ, against a CMN 
paramagnetic salt in zero field for temperatures below 0.1K and down to 7 mK, and against 
reference thermometers (calibrated with the CMN) placed in the compensated zone of the 
magnet for field measurements. The electrical resistivity was measured simultaneously, with 
a conventional lock-in detection at low frequency (around 2 Hz, for minimum quadrature 
signal), allowing a precise determination of L(T). 

Figure [l] compares the electrical resistivity p, and the thermal resistivity w t h = 
at low temperature (below 0.15K), for zero magnetic field (H = OT) and very slightly 
above the critical field H c (H = 70mT). This is the most important result, but let us 
first describe the insert which gives a broader view thanks to an extended temperature 
range (from 7mK up to 7K). It displays the Lorenz ratio L(T)/Lq for both fields, with a 
logarithmic temperature scale: There is a large temperature range, from 3 K down to 0.02 K, 
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where the normalized ratio L(T)/Lq is clearly below 1, down to a value L(T)/Lq fa 0.87. 
This is indicative of dominant inelastic scattering (the so-called “vertical processes”), which 
affect more strongly thermal than electrical transport. The rather unusual feature is that 
despite a large residual resistivity (the RRR is “only” 65 for our sample), inelastic scattering 
has a significant contribution down to such low temperatures. But this is consistent with the 
marked temperature dependence of the resistivity down to the lowest measured temperatures 
(8mK), as can be seen on the main graph. Inversely, the fast increase of the Lorenz ratio 
above 2 K, overpassing 1 above 3 K, indicates additional contributions to the heat transport 
(like phonons). 

These “bosonic” contributions are expected to vanish rapidly below 1 K, and therefore, 
the new increase of the Lorenz ratio above 1 below 20 mK is again a very unusual feature, 
which could not be anticipated from the resistivity behavior. This low-temperature increase 
of the Lorenz ratio has been already observed below 30 mK by Pfau et al. [5], on a slightly 
more resistive sample (RRR fa 45 instead of 65). But we insist that this is probably the most 
puzzling feature of thermal transport in YbRl^S^: How is it that new bosonic excitations, 
able to contribute as a new heat channel, appear at such low temperatures? Moreover, 
contrary to the statement in [5], we observe that this increase is essentially unchanged for a 
field of 70 mT, suppressing the long range AF order, and the putative magnon (long lived) 
excitations expected in this ordered phase. Discussion of the origin and field evolution of 
this contribution is continued later in this paper. 

The main plot of Fig. [l] shows the detailed evolution of p and w t h below 0.15 K, down 
to 10 mK: In particular, one clearly observes a kink at T/v in zero held on the electrical 
and thermal transport. It is suppressed by an applied held of 70 mT (very close to H c ), 
which also lowers both resistivities. It demonstrates that the critical fluctuations emerging 
in the neighborhood of the AF transition contribute to the scattering of the quasiparticles. 
However, this scattering seems to be rapidly damped below T N , as no change of w t h is 
observed below 50 mK, between zero held and 70mT. Concerning now the additional low 
temperature contribution, the strong decrease of w t h below 30 mK appears independent of 
the applied held. It already shows the difficulty to discuss the validity of the WFL: The linear 
extrapolation of the data of w t h between 30 and 50 mK is compatible, within the dispersion 
of the points, with that of the resistivity at OK. The whole question is to control these 
extrapolations, or to estimate the amplitude and range of the extra contribution responsible 
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for the low temperature drop of w t h • But in any case, on the bare data, no drastic change 
is seen between zero held and H c . 

Figure [2] compares the relative temperature variation of L(T)/L 0 for holds above H C1 
up to 4T (see Supplementary Material A, figure [5] for the data of k, and p). It appears 
that the low temperature increase of the Lorenz ratio is held sensitive above 100 mT, and 
rapidly suppressed for rather weak helds (above 0.2 T). This supports a magnetic origin for 
this additional heat channel, which is reasonable as we did not expect a sizable phonon 
contribution at such low temperatures. In parallel, we do observe that the decrease below 
one of the Lorenz ratio above 30 mK is also progressively suppressed under held, pointing 
to a magnetic origin of the low temperature inelastic scattering. However under held, for 
H > 0.5T (see also insert up to 4T), L(T —* 0) seems to saturate at 0.97 L 0 . This small 
deviation from may be due to experimental precision, or related to the complex multiband 
structure of the FS. In any case, the fact that this value does not change with helds much 
larger than H c means that it cannot be seen as a violation of the WFL at the QCP. In the 
following, particularly to estimate the electronic contribution, the Lorenz number is taken 
at this value found for T —> OK: L(0) = 0.97Lo- 

As regards the question of the origin of the additional contribution to thermal transport 
appearing below 30 mK (hereafter labelled K a dd ), all hints point to a magnetic origin: (i) 
even if it survives the critical held, it is rapidly suppressed by helds of order of 0.5 T; (ii) 
a new phononic contribution is completely unexpected in this temperature range: It would 
require a strong lattice instability, which has not been detected so far; (iii) besides the 
antiferromagnetic long range order at 70 mK, other magnetic transitions have been reported 
at 12mK |15j and/or 2mK [16]. In reference [5], n a dd has been attributed to magnons, 
an assumption supported by comparison to the specihc heat ( C p ) data: The temperature 
dependence of C p has a bare T 3 regime below T* = 55 mK [3J [5], which is attributed to 
AF magnons. Because it exists even when the long range ordered state is suppressed, we 
would rather compare n a dd in YbRl^SL to a similar additional contribution to heat transport 
observed in the ferromagnet UCoGe, both above and below T Curie, and also in a very low 
temperature range (for the heat current along the b-axis) [IT]. Obviously, “paramagnons” or 
overdamped spin huctuations, can also contribute to an additional heat transport channel. 

Quantitatively, we will perform a similar analysis as in HU, assuming that we can separate 
the electronic contribution in a quasiparticle K qv channel, which should follow (or not!) the 
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FIG. 2. (color online) Ratio L(T)/Lq at low temperatures below IT. The dashed lines are the 
extrapolation to L(0) = 0.97Lo when T > 30mK. In inset, L(T)/Lq at higher fields, which extrap¬ 
olates also to 0.97 at low temperature. 



FIG. 3. (color online) Comparison of the electronic thermal conductivity using the expression [T] 
with a = 1 and a = 2.2. The latter value was chosen to optimize the low temperature contribution 
of Kqp. 
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WFL, and a “spin-fluctuation” contribution n a dd , expected for example to be much less 
sensitive to the RRR. As for the case of UCoGe, a difficulty for quantitative estimation of 
K a dd is that the quasiparticle contribution is badly known, due to the importance of inelastic 
scattering even at 30 mK (see the Lorenz ratio in figure [TJ. So we assume that n qp can be 
described by [18] : 

K qp _ _ L(0) _ j 

T p 0 + a(p — p 0 ) 

This expression takes into account the stronger effect of the electronic inelastic scattering 
on the thermal transport (than on the electrical transport), through the factor a. po is 
the residual resistivity and a = 1 + W vert /W} lor , where W vert and Wh or are respectively the 
scattering rates due to vertical (small wave-vector transfer q, inelastic) and horizontal (large 
q, mostly elastic) processes, assuming that Mathiessen’s rule holds. The crude simplification 
is to take a constant, which can be valid for temperatures much smaller than the typical 
energy of the fluctuations responsible for the inelastic scattering. This simplification has 
already been widely used for the heavy fermion systems like UPt 3 [T9], CeRhIn 5 [20] and 
UCoGe [T7]. The figure [3] compares the bare data with different estimates of K qp for a = 1 
i.e L(T ) = L(0) and a = 2.2, at 0T (in the AF state), 70mT (close to H c ) and 0.27T (in 
the PM state). Clearly, a = 1 overestimates the quasiparticle contribution, whereas a = 2.2 
(the best value we found) gives an estimate of K qp compatible with the data for the three 
fields. It even gives a fair account of the data between 30 mK and 100 mK at 0.27 T. The 
extracted extra contribution is shown in the Supplementary Material A, figure [6j 

The most important information from the analysis shown on figure [3] is that little change 
is observed between zero field and H c , apart from the suppression of the anomaly at T^. 
Otherwise, (i) it is possible to analyze the data of thermal transport in YbRh 2 Si 2 , assuming 
that the quasiparticle contribution follows the WFL even at H c , and with the same coefficient 
a = 2.2 describing the effect of inelastic scattering; (ii) K a dd is little influenced by the 
suppression of the Neel temperature at H c , and really starts to decrease only above this field 
(see also figure [6] of the Supplementary Material A), pointing to a contribution of driven 
(robust) short range correlations. 

For point (i), it will be very difficult to reach firmer conclusions, as even lower temperature 
measurements would face the problem of the dominant additional contributions, and of the 
difficulty of a robust quantitative analysis, together with the regimes introduced by the 
additional transitions at 12 and 2 mK nans]. For point (ii), this is similar to observations 



on the thermopower, where the change of regime under field occurs far above H c M, and 
little is seen when crossing H c . 

Moreover, the question of the role of magnetic fluctuations in the heat transport in 
YbRli 2 Si 2 is not limited to the low temperature regime. Also at high temperatures, above 
IK, we can detect a field dependence of w t h , which (like for resistivity), points to a contri¬ 
bution of magnetic fluctuations on the scattering mechanisms. In this temperature range, 
microscopic measurements have revealed a complex situation: NMR experiments [2JLJ have 
shown that AF fluctuations with finite wavevectors compete with ferromagnetic (FM) fluc¬ 
tuations; recent investigations by NMR [Tj gives the image of interfering Fermi liquid and 
non Fermi liquid component close to H c . Neutron scattering experiments [2] down to 0.1 K 
up to 10 T along the c-axis confirm that incommensurate AF correlations become dominant 
on cooling against the FM correlations, while the magnetic field induces a Zeeman resonant 
excitation. These FM fluctuations could be an important ingredient to understand the ob¬ 
servation of a well defined electron spin resonance in YbRl^SL [23, j25|. So at least far from 
T/v or at high magnetic field, FM coupling may be dominant, pointing to the proximity of a 
FM QCP [26], 27]. We do find (see Supplementary material B) that the thermal resistivity 
above IK behaves as predicted by Ueda and Moriya [3] when inelastic scattering is controlled 
by ferromagnetic fluctuations. 

To summarize we have investigated in detail the thermal transport of YbRh 2 Si 2 , in the 
AF phase and in the PM phase, below and above H C1 and extended the measurements down 
to 10 mK. This study clearly shows that an additional magnetic contribution to thermal 
conductivity appears below 30 mK, with a gradual decrease for fields above the critical field 
H c . This has been deduced with an electronic quasiparticle contribution, which also has no 
peculiar singularity at H c , and satisfies the WFL even at H c . To conclude unambiguously, 
by direct measurements, if the WFL is satisfied or violated, would require measurements 
down to the mK temperature range, where the additional contribution should eventually 
disappear. But other recent experiments also challenge the “local QCP scenario”, which 
could have led to the violation of the WFL at the critical field: For example, ARPES 
measurements [29., 30] show that above T/v, the 4f-clectrons in the paramagnetic phase are 
clearly itinerant, which leaves open the validity of the image of small and large Fermi surface 
fluctuations developed in the Kondo breakdown scenario. A new appealing approach is the 
so-called “strong coupling theory” of heavy-fermion quantum criticality, where large critical 


9 


spin fluctuations at the AF wave-vector induce fluctuations at small wave-vectors, producing 
a diverging effective mass over the entire Fermi surface [31]. It would be interesting to see 
if these fluctuations could be responsible for the detected low temperature contribution, as 
well as to the large inelastic scattering observed down to very low temperatures ? 

This work has been supported by the French ANR grant SINUS and the ERC grant 
“NewHeavy Fermion”. 
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Supplementary Materials 


SUPPLEMENTARY MATERIAL A: DETAILS ON THE TRANSPORT MEASURE¬ 
MENTS 

As a complement to figure [T| of the main paper, displaying the thermal and electrical 
resistivities below 0.15 K, figure [4] presents p and w t h at 0 T up to 7 K: This shows the strong 
temperature dependence of p in the whole temperature range, as well as the saturation of 
Wth above 3K, due to the growing contribution of bosonic excitations to heat transport. 



FIG. 4. (color online) Electrical and thermal resistivities at 0T in the whole temperature range 
from 10 mK up to 7 K. 

If several (parallel) channels contribute to the thermal transport, it is better to discuss 
the thermal conductivity rather than w t h■ It is displayed on figure [5j at low temperature 
and up to 0.1 T, and the corresponding electrical resistivity is shown in the inset. The 
smooth variation of the thermal conductivity through H c is confirmed by the temperature 
dependence of n/T below 100 mK on each side of H c for different magnetic fields. 

As explained in the main paper, we could analyze the quasiparticle contribution, tak¬ 
ing into account the peculiar effects of inelastic scattering on heat transport, and with no 
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FIG. 5. (color online) Thermal conductivity of YbRh 2 Si 2 at low temperature below 0.1 T. In inset, 
the corresponding electrical resistivity. The dash line shows the Neel temperature at zero field. 

detected anomaly when crossing the critical field H c . Figure [6] presents the additional con¬ 
tribution Kadd/T for three different fields H = 0, H ~ H c and H ~ 4 H c , subtracting this 
calculated value of n qp (T). As announced, there is little change of K a dd at H c , whereas its 
collapse is noticeable for a field of 0.27T. For even higher field values, K a dd rapidly falls 
bellow the error bars of the analysis, an effect reinforced by the increase of the bottom 
temperature of the measurements due to prohibitive equilibrium relaxation times. But as 
regards the direct test of the validity of the WFL, the strong “divergence” of K a dd/T at low 
field at 10 mK sets the level of the experimental challenge: Probably temperatures of order 
a few mK at most, are required to observe the downturn of this contribution? 


SUPPLEMENTARY MATERIAL B: DETAILS ON THE HIGH TEMPERATURE 
ANALYSIS 

As underlined, evidences exist that the main magnetic intersite interaction is tempera¬ 
ture dependent. As the “average” effective mass increases drastically on cooling, the rela¬ 
tive weight of the different outbands may change and produce a feedback on the magnetic 
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FIG. 6. (color online) Estimation of the very low temperature extra contribution to the thermal 
conductivity at OT, 70mT and 0.27 T. No significant changes are seen between 0T and 70 rnT, 
indicating that n a dd has certainly the same origin at both fields. n a dd exists for fields higher than 
H Cl but with a decreasing amplitude. Lines are guides to the eyes. 

coupling itself. Due to evidences for FM fluctuations above 1K DEI, we investigated if 
YbRh 2 Si 2 could be seen, far above T N , as a weak ferromagnet above its Curie temperature 
from the point of view of thermal transport? In this case, if only the quasiparticle-magnetic 
fluctuations interactions are taken into account, Ueda and Moriya have predicted that the 
thermal resistance should behave as [3]: 


wth = P + BT (2) 

So, even if the electrical resistivity has no real simple temperature dependence (see figures 
3-4 of reference [3], the additional thermal resistance for weak ferromagnets above Tc U rie 
or near ferromagnets at high enough temperatures is roughly linear in temperature (with 
our definition of w t h)'- See figure 5 of [3]. Naturally, in this temperature range, we also 
expect a phononic contribution, limited by electron-phonon scattering, which should follow 
a quadratic temperature dependence: k p h = PT 2 , with P temperature and field independent. 
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We obtain thus for the thermal conductivity at high temperature: 



( 3 ) 

( 4 ) 


Once the phonons are subtracted, the thermal resistance w t h = u) qp = represents the 


electronic heat channel with the scattering of quasiparticles by FM fluctuations. We call 
the additional linear term in the thermal resistivity w qp _fi uc = BT. Contrary to the low 
temperature situation, the magnetic fluctuations above 1K do not contribute with their own 
channel to heat transport: They only have the “usual” contribution to the inelastic scattering 
of the quasiparticles. The expression [4] is very simple, as p is determined experimentally and 
the parameter P is held independent. Only the parameter B can vary with the held. Such a 
decomposition of the thermal conductivity has already been successfully used for the weak 
ferromagnet ZrZii 2 |4j. 

Surprisingly, after removal of the phonon and resistivity contribution, we do find that 
Wqp-fiuc is linear above 1 K, for helds up to 4T, as shown in figure [7J for a held independent 
parameter P = 0.022W.K _3 .m _1 , close to 0.017W.K _3 .m _1 used by Pfau et al. |5|. At 
lower temperatures, a deviation of w qp _fi uc from linearity appears, which may be due to the 
growing importance of AF fluctuations. 

The evolution of the B coefficient is shown in figure [8j It increases linearly with the mag¬ 
netic held, as if the applied held and the warming reinforce the ferromagnetic huctuations. 
This behaviour is quite singular and may be related to the intermediate valent character of 
YbRh 2 Si 2 . In Ce compounds, like in the CeRu 2 Si 2 family, ferromagnetic huctuations are 
just enhanced at the metamagnetic held where simultaneously the AF interactions drops [6]. 


* Present address: Low Temperature Laboratory, Aalto University, P.O. Box 13500, FI-00076 
Aalto, Finland 
1 jean-pascal.brison@cea.fr 
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FIG. 7. (color online) Electronic thermal resistivity w qp _fi uc up to 5K after subtraction of the 
electrical resistivity. For clarity, an offset has been added to each curve (+0.2^12.cm at each held 
step). Dashed lines are guides to the eyes. 
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FIG. 8. (color online) Evolution of the coefficient B with the magnetic held. 
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